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Abstract 

Critical properties of the compact three-dimensional U(l) lattice gauge theory 
are explored at finite temperatures on an asymmetric lattice. For vanishing value 
of the spatial gauge coupling one obtains an effective two-dimensional spin model 
which describes the interaction between Polyakov loops. We study numerically the 
effective spin model for Nt = 1,4,8 on lattices with spatial extension ranging from 
L = 64 to L = 256. Our results indicate that the finite-temperature U(l) lattice 
gauge theory belongs to the universality class of the two-dimensional XY model, 
thus supporting the Svetitsky-Yaffe conjecture. 
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1 Introduction 



The finite temperature behavior of the compact three-dimensional (3d) U(l) lattice gauge 
theory (LGT) is the subject of numerous investigations (see, e.g., Ref. pQ and references 
therein). It is well-known that at zero temperature the theory is confining at all values of 
the bare coupling constant [2]. At finite temperature the theory undergoes a deconfine- 
ment phase transition. Both phenomena are expected to take place in Ad QCD as well. 
Therefore, the 3d U(l) gauge theory constitutes one of the simplest models with contin- 
uous gauge symmetry which possess the same fundamental properties as QCD. In view 
of these common features the critical properties of 3d U(l) LGT deserve comprehensive 
qualitative and quantitative understanding. 

On the theoretical side one should mention two results regarding the critical behavior 
of 3d U(l) LGT. The first result states that the partition function of 3d U(l) LGT in 
the Villain formulation coincides with that of the 2d XY model in the leading order of 
the high-temperature expansion [3]. In particular, the monopoles of the original U(l) 
gauge theory are reduced to vortices of the 2d system. The second result follows from 
the Svetitsky-Yaffe conjecture: the finite-temperature phase transition in the 3d U(l) 
LGT should belong to the universality class of the 2d XY model if correlation length 
diverges [3|. Then, two possibilities arise: either the transition is first order or it is the 
same transition which occurs in the 2d XY model. The XY model is known to have 
Berezinskii-Kosterlitz-Thouless (BKT) phase transition of infinite order [51 [5]. Several 
important facts could be deduced from these results. First of all, the global U(l) symmetry 
cannot be broken spontaneously even at high temperatures because of the Mermin- Wagner 
theorem. Consequently, a local order parameter does not exist. Secondly, one might 
expect the critical behavior of the Polyakov loop correlation function T(R) to be governed 
by the following expressions 

r(«) x , (i) 

for (3 > (3 C and 

T(R) x exp [-R/m\ , (2) 

for (3 < /3 C , t = (3 C / (3 — 1. Here, R ^> 1 is the distance between test charges and £ ~ e bt 
is the correlation length. Such behavior of £ defines the so-called essential scaling. The 
critical indices n(T) and v are known from the renormalization-group analysis of the XY 
model: i](T c ) = 1/4 and v = 1/2, where T c is the BKT critical point. Therefore, the 
critical indices n and v should be the same in the finite-temperature U(l) model if the 
Svetitsky-Yaffe conjecture holds in this case. 

The first renormalization-group calculations of the critical indices, presented in [3], 
gave support to the conjecture even though they did not constitute a rigorous proof. 
The direct numerical check of these predictions was performed on lattices N% x N t with 
N s = 16, 32 and Nt = 4, 6, 8 in Ref. [?]. Though authors of Ref. [7j confirm the expected 
BKT nature of the phase transition, the reported critical index is almost three times that 
predicted for the XY model, n(T c ) w 0.78. More recent numerical simulations of Ref. [1] 
have been mostly concentrated on the study of the properties of the high-temperature 
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phase. We have to conclude that, so far, there are no numerical indications that the 
critical indices of 3d U(l) LGT do coincide with those of the 2d XY model. Moreover, 
since a rigorous determination of the critical indices is not available even for the XY 
model one can hardly hope for a rigorous analysis of the critical behavior of 3d U(l) 
LGT. 

The absence of reliable results in the vicinity of the BKT critical point was our primary 
motivation to study the deconfmement phase transition in 3d U(l) LGT. The difficulties 
in computations of critical indices of the XY model are well-known and we do not intend 
to discuss them here (see Ref. [8] for a summary of recent results and problems). It should 
be clear however, that in the context of the 3d theory a reliable determination of critical 
properties becomes even harder and requires simulations on very large lattices. We have 
decided therefore to attack the problem in a few steps. Consider the finite-temperature 
model on anisotropic lattice with different spatial and temporal coupling constants; as 
a first step, in this paper we investigate the limit of vanishing spatial coupling. The 
major advantage of this limit is that the integration over spatial links can be performed 
analytically. The result of such integration is an effective two-dimensional spin model for 
the Polyakov loops. The latter can be studied numerically. 

This paper is organized as follows. In the next section we introduce the compact U(l) 
LGT on anisotropic lattice and study it for vanishing spatial coupling. In the Section 3 
we describe briefly our numerical procedure. The result of simulations are presented in 
the Section 4. Conclusions and perspectives are given in the Section 5. 



2 The 3d U(l) lattice gauge theory 

We work on a 3d lattice A = L 2 x N t with spatial extension L and temporal extension N t . 
Periodic boundary conditions on gauge fields are imposed in all directions. We introduce 
anisotropic dimensionless couplings in a standard way as 

A = -5— , Ps = -4- = A e , £ = - , (3) 

g a t g A a s a s 

where a t (a s ) is lattice spacing in the time (space) direction, g 2 is the continuum coupling 
constant with dimension a -1 . The finite-temperature limit is constructed as 

£ _ o , N t , L -> oo , a t N t = ^, (4) 

where T is the temperature. 

The 3d U(l) LGT on the anisotropic lattice is defined through its partition function 

as 



where S is the Wilson action 

S[u) = P„y] cosuj(p s ) + Pt^2 COSUJ (Pt) ( 6 ) 



pt 



3 



and sums run over all space-like (p s ) and time-like (p t ) plaquettes. The plaquette angles 
uj{p) are denned in the standard way. The correlation of two Polyakov loops can be 
written as, e.g. 



T(R) 



cxp 



N t -1 

E M)Oo,£i- 





X2) 



u (x ,xi,x 2 + R)) 



.to 



(7) 



As stated in the Introduction we would like to explore the limit fi s = 0. Consider the 
strong coupling expansion at fi s <C 1. The general form of such expansion reads 



Z(fi u fi s ) = Z{fi u fi s = 0) Z 2k (fi t 



k=l 

In this paper we study the zero-order partition function Z(fi t , fi s = 0) defined below. The 
series on the right-hand side of the last expression is known to be convergent uniformly 
in the volume, both for the free energy and for the gauge-invariant correlation functions. 
The uniform convergence guarantees the existence of the limit N t —>■ oo. The strong 
coupling expansion, done even in one parameter, might be far from the continuum limit. 
Nevertheless, one expects that already the zero-order approximation captures correctly 
the critical behavior of the full theory. An example is given by the following Polyakov 
loop model 

S eS = PeffJ2 Re W(x)W*(x + n) , (9) 



derived at finite temperature for (d + l)-dimensional SU(N) pure gauge theory in the 
limit fi s = 0. Here, (3 e g oc (3 t *. As is well known, this model reveals correctly the critical 
behavior of the original theory thus supporting our approximation. 

In the zero-order approximation the integration over spatial gauge fields can be easily 
done and leads to the following expression for the partition function 



z(p t ,p s = o) 



2tt 



E 



I?* (Pi) exp{ir(u(x)-uj(x + e n ))} 



(10) 



where x belongs to the two-dimensional lattice A 2 = L 2 and u>(x) = u(xi, x 2 ). Here, I r (x) 
are modified Bessel functions and e tTU) ^ is the Polyakov loop in the representation r. 
For N t = 1 using the formula ^ r I r (x)e truJ = e XCOSUJ one finds 



= 0)| 



N t =l 



2tt 

n 



dw(x) 
2rt 



cxp 



A E cos^^) - u(x + e n )) 



which is the partition function of the 2d XY model. Thus, in this case the dynamics of 
the system is governed by the XY model with the inverse temperature fit. For N t > 2 
the model (FlOi) is of the XY-type, i.e. it describes interaction between nearest neighbors 
spins (Polyakov loops) and possesses the global U(l) symmetry. Moreover, consider now 
two different limits - the strong coupling limit fi t <C 1 and the weak coupling limit fi t ^> 1. 
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In the leading order of the strong coupling limit one can easily find from (fIU|) . up to 
an irrelevant constant, 



z(A«i,ft = o) 



2vr 

n 



duj(x) 
2tt 



cxp 



h(f3 t ) cos(cj(:e) — u(x + e n )) 



'121 



which is again the XY model with the coupling h given by 



The Polyakov loop vanishes while the correlations of the Polyakov loops are given, at the 
leading order, by 



r(R) 



R 



(13) 



To study the weak coupling limit it is convenient to perform duality transformations 
which are well-known for the XY model. Taking then the asymptotics of the Bessel 
functions one obtains, up to an irrelevant constant, 



Z(A>1,& = 0) 



ex p 

r(x)=— oo 



1 s 



x n=l 



(14) 



This is nothing but the Villain version of the XY model in the dual formulation with an 
effective coupling 

P = N t /p t = g 2 /T . (15) 



This shows that the region (3 a = 0, (3 t 3> 1 is also described by the XY model. 
In the general case of arbitrary (3 t the full effective action 



S e fi = CkCosk(u(x) -uj(x + e n )) 



(16) 



x,n k 



will include all representations k of the Polyakov loops. In our case the coefficients Ck are 
given by 



i>2-k J °° 

C k = / — cosku log{l + 2V[6 r (A)]^cosrw) 
Jo 2tt I ^ J 



(17) 



where b r (Pt) = Ir{flt) / lo(Pt) ■ If there is a critical point at which the correlation length is 
divergent then on general grounds (universality, limiting behavior, etc.) one assumes that 
the model described by the effective action (flBTl indeed possesses the same critical behavior 
as the XY model. Nevertheless, we are not aware of any direct numerical check of the 
universality for models of the type (fl6|) if C k ^ for k = 2,3, .... In the following sections 
we present numerical simulations which give support for the expected BKT behavior of 
the model (fTBTl . Our results hold only for the model with Ck defined by (TT7I) . We would 



5 



like to stress that it is not obvious that for all possible the correlation length really 
diverges. For example, it was proven in Ref. [H] that the model with coefficients 

with sufficiently large p, exhibits a first order phase transition, so that one could expect 
that the correlation length stays finite across the phase transition point. 



3 Numerical set-up 

Determining the universality class of the 3d U(l) gauge theory discretized on a L 2 x N t 
lattice means determining its critical indices. A convenient way to accomplish this task 
is to study the scaling with the spatial size L of the vacuum expectation value of suitable 
observables, determined through numerical Monte Carlo simulations. 

For the special case (3 a = 0, one can take advantage of Eq. (ITOl and describe the original 
gauge system with a two-dimensional spin model whose action 5" is defined through 

z(A,/?. = o)= riIi7 ex P 5 ' ( 19 ) 

" X 

and reads 

oo 

S' = J2 iogj 1 + 2 cosr(uj(x) - u(x + e n ))} . (20) 

x,n r=l 

The infinite series in r can be truncated early, since the 6 r 's vanish very rapidly for 
increasing r. We studied the dimensionally reduced system with the Metropolis algorithm, 
taking the first twenty b r couplings (notice that 6 2 o (A = 1) ~ 1CT 25 ). 

Our goal is to bring evidence that the system exhibits BKT critical behavior for any 
fixed N t . This is trivially verified in the case N t — 1, since by inspection of Eqs. ffT91 
and (|2"0"|) . the theory reduces exactly to the XY model. Therefore the case N t = 1 can be 
used as a test-field for the description and the validation of our procedure. 

Before presenting numerical results it is instructive to give some simple analytical 
predictions for the critical values (3 t at different values of Nt- Such critical values can be 
easily estimated if one knows /3^ r for N t = 1. Since the model with Nt — 1 coincides with 
the XY model one has j3^ r (N t = 1) ~ 1.119 and approximate critical points for other 
values of N t can be computed from the equality 

&i(l.H9) = [bM r )] Nt ■ (21) 

Solving the last equation numerically one finds /5 t cr . The results are given in the Ta- 
ble [TJ As will be seen below the predicted values are in a reasonable agreement with the 
numerical results. 



6 



Table 1: Analytical estimates of /5 4 cr for several values of N t (first row) compared with the 
numerical results obtained in Section H] (second row). 





2 


4 


8 


16 




2.0003 


3.39389 
3.42(1) 


6.10642 
6.38(5) 


11.6385 



4 Results at /3 S = 

4.1 N t = 1 

The main indication of BKT critical behavior is a peculiar scaling of the pseudo-critical 
coupling with the spatial lattice size L, consequence of the essential scaling, 

^-A-OSSIF" (22) 

where (3 pc (L) is the pseudo-critical coupling on a lattice with spatial extent L, (3 C is the 
(non-universal) infinite volume critical coupling and v is the (universal) thermal critical 
index. 

The pseudo-critical coupling f3 pc (L) is determined by the value of j3 for which a peak 
shows up in the susceptibility of the Polyakov loop, 

X = L 2 (\P\ 2 } 1 ^ = 7^E P - ( 23 ) 

X 

here the local Polyakov loop variable P x corresponds to the spin s x = expiu(x) of the 
XY model. In Fig. [1] we show the behavior of the absolute value of the Polyakov loop \P\ 
(top) and of the susceptibility x (bottom), for varying (3 on lattices with L— 32, 64, 128. 

To extract (3 pc (L) in a more reliable way, we performed the multi-histogram inter- 
polation [10]; errors were determined by the jackknife method. Results for /3 pc (L) are 
summarized in Table [2j 

Table 2: (3 pc {L) for N t =l, 4, 8 and for several values of L. Errors are determined by a 
jackknife analysis. 



L 


N t = l 


N t = A 


N t = 8 


64 




3.1250(51) 


5.531(19) 


128 


1.0051(16) 




5.754(22) 


150 


1.0094(26) 


3.2190(40) 


5.7945(59) 


200 


1.0227(15) 


3.2368(39) 




256 


1.0278(20) 







1 Throughout this Section we use the notation (3 t = /3. 
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We determined (3 c (N t = 1) by fitting the pseudo-critical coupling /3 pc (L) given in the 
second column of Table [2] with the law 

« L » = A+ (bi7)W < 24 > 

in which v was fixed by hand at the XY value, v = 1/2. We got fl c (N t = 1) = 1.107(9) 
and A(N t = 1) = —2.4(2) (x 2 /d.o.f.=0.78), which is quite in agreement with the best 
known XY critical coupling, (3 C = 1.1199(1), given in Ref. [TT] . 

The determination of (3 C is crucial in order to extract critical indices; indeed, they 
enter scaling laws which hold just at /3 C , such as, for example, 

X m-L 2 -^ , (25) 

where 7] c is the magnetic critical index. Actually in Eq. [25] one should consider logarithmic 
corrections (see [121 EE] and references therein) and, indeed, recent works on the XY 
universality class generally include them. However, taking these corrections into account 
for extracting critical indices calls for very large lattices even in the XY model; for the 
theory under consideration to be computationally tractable, we have no choice but to 
neglect logarithmic corrections. 

We determined x(P — 1-12) for L=64, 128, 150, 200, 256 - see Table [3]for a summary 
of the results. Fitting with the law f l25l) . we found r\ c = 0.256(29) (x 2 /d.o.f.=0.2), in nice 
agreement with the XY value, t] c = 1/4. The same analysis repeated at (3 = 1.107, i.e. 
at the central value of our determination of /3 C , on lattices with L=64, 128, 200, gave 
r] c = 0.237(61) (x 2 /d.o.f.=0.01). 

An alternative strategy to determine i] c uses the large distance behavior of the point- 
point correlator of the Polyakov loop, 

C(R) = J2^(PlP*+Re n ) , (26) 

x,n 

where e n is the unit vector in the n-th direction. Without logarithmic corrections, one 
has 

C(R)~^r c ■ (2T) 

In Fig. [2] we plot log C(R) vs \ogR for L=200 at (3 = 1.12; linearity is clear up to 
R ~ 30. Deviations at larger distances are due to finite size effects (echo terms are 
expected to be strong, since the correlator is long-ranged) and possibly to logarithmic 
corrections. In the linear regime (5 < R < 30), the naive fit with a power law gives 
7]= 0.22942(31) (x 2 /d.o.f.=0.83). The same analysis at (3 = 1.107 and L=200 gives 7]= 
0.2380(20) (x 2 /d.o.f.=0.05) in the range 1 < R < 45. On the same volume one sees that, 
for lower /3's, rj goes towards the expected value and that the linear region gets wider and 
wider. 

The effective rj c index, defined as 



s 



must exhibit a plateau in the region where (1271) holds. Fig. E^top) shows that the larger 
the volume the larger the region in which there is a plateau at small distances. The chosen 
value of Rq must belong to the linear region in order to minimize finite size effects. We 
have verified that varying Rq in the linear region does not change the result and have 
chosen i?o = 10 for all the cases considered here. 

Since for the larger lattices (L=200 and L=256) plateaux are overlapping at small 
distances, one can conclude that thermodynamic limit is reached. We estimate the plateau 
value from the most precise data we have (L=200) as r]((3 = 1.12) = r] c s(R = 6) = 
0.23101(49), since the latter is the value of r) e s in the linear region compatible with the 
largest number of subsequent points. Deviations from the expected value ^=0.25 can be 
due either to logarithmic corrections or to the overestimation of (3 C . Repeating the same 
procedure for slightly lower (3's we find: r]((3 = 1.115) = 0.23491(47) and r]((3 = 1.107) = 
0.24085(44). Notice that rj approaches the expected value as (3 lowers. The relation 
between r\ and (3 is well described by a linear function (x 2 /d.o.f.=0.04) and this suggests 
that the (3 value at which 7?=0.25 is really close to those considered. Fig. [3](bottom) shows 
the correlation function C(R) rescaled by L~ v in units of R/L; it turns out that, when 
the best determination for rj is used (in the present case, rj = 0.23101) data from different 
lattices fall on top of each other over almost all the range of distances considered. 

There are other observables which turned out to be useful in establishing the BKT 
scaling in the 2d XY model and which we do not use in the present work: the helicity 
modulus y [U1[T3], the second moment correlation length £ 2 (see, for instance, |13j) and 
the U4 cumulant, proposed in [IB] . We plan to use them all when we will study the 
general case (3 S 7^ 0. For the purposes of the present work we have only tried to use the 
U4 cumulant, but both lattice sizes and statistics seem to be not enough large to extract 
any useful information from this observable. 

4.2 N t =4 and 8 

In this Subsection we extend the study performed in the N t = 1 case to the cases of N t = 4 
and 8, with the aim of showing that the universal XY features are not lost increasing N t 
at p s =0. 

In Table [2] we give the values of the pseudo-critical couplings f3 pc (L) obtained from the 
peaks of the Polyakov loop susceptibility for several values of L at N t = 4 and 8. Fitting 
these values with the law (124"]) with v = 1/2 fixed, we get 

(3 c {N t = A) = 3.42(1), A(N t = 4) = -5.1(3) , ( X 2 /d.o.f. = 0.43) 

(3 c (N t = 8) = 6.38(5), A(N t = 8) = -15(1) , ( X 2 /d.o.L = 0.006) . (29) 

This result shows that essential scaling is satisfied, i. e. in both cases transition is compat- 
ible with BKT. It is worth noting that these values of j3 c are in nice agreement with the 
estimates given in Table [lj This suggests that the dynamics of the effective model near 
the transition point is indeed dominated by the lower representations, thus justifying the 
truncation of the series in Eq. (I20j) . 
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Table 3: x{L) for N t =l, 4, 8. Errors are determined by a jackknife analysis. 



L 


N t = l 


N t = A 


N t = 8 


64 


7.19(12) 


9.30(57) 


7.33(37) 


128 


24.6(1.3) 


35.9(4.4) 


25.1(1.5) 


150 


32.9(1.9) 


42.5(2.5) 


32.5(1.9) 


200 


51.4(2.7) 


65.2(2.8) 


58.4(3.3) 


256 


80.3(4.2) 


101.7(5.4) 


86.4(3.6) 



In Table |3] we give the values of the Polyakov loop susceptibility for several values of 
L at (3 = 3.42 for N t = 4 and at (3 = 6.38 for N t = 8. Fitting with ([25]), we find 

Vc (N t = 4) = 0.290(54) (x 2 / d - oi - = °- 69 ) 

rj e (N t = 8) = 0.212(46) (x 2 /d.o.f. = 0.28) . (30) 

Results agree with the universal XY value r\ c — 1/4, although errors are quite large. 

A more precise determination of the magnetic index can be achieved through the 
study of the point-point correlation function. In Figs. BJtop) and[5](top) we show r] c s(R) 
for three values of the spatial size L for the cases of N t = 4 and N t = 8, respectively. Our 
estimated plateau values, taken from data at L = 200, are 

7^/3 = 3.42) = V es(N t = 4, R = 2) = 0.2724(11), 
T)((3 = 6.38) = rj eS (N t = 8, R = 3) = 0.2499(11) . 

For N t = 4, rj overshoots by little the XY universal value, while for Nt — 8 it is in 
nice accord with it. The deviation for Nt — 4 is most likely washed out by a fine tuning 
of the critical coupling within its error bars. 

One can observe, moreover, that the shape of the curve of values of r) e g(R) changes 
qualitatively in the same way when the thermodynamic limit is approached for N t = 1 
and Nt = 8, while it has a different behavior for N t = 4. This may be an indication that 
for Nt = 1 and Nt = 8 at the /3's chosen for the simulation the system is in the same 
phase (P > (3c), i.e. correlators have the same behavior. 

Figs. H^bottom) and [^bottom) show the correlation function C(R) rescaled by L~ n 
in units of R/L, with r\ fixed at the central value of our determinations (77 = 0.2724 for 
Nt = 4 and rj = 0.2499 for N t = 8); one can see that data from different lattices fall on 
top of each other over a wide range of distances. 

In summary, essential scaling is verified both for N t = 4 and 8, thus indicating that 
indeed the occurring transitions are compatible with BKT. Moreover data point to values 
of the thermal and magnetic critical indices of the 2d XY universality class. This leads 
us to conclude that for N t = 4 and 8 the 3d U(l) LGT at f3 s =0 belongs to the 2d XY 
universality class and this supports the conjecture that the same holds, in general, for 
any N t at (3 S = 0. 
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Since we do not study the correlation length, we are not allowed to rule out the 
possibility that it keeps finite and the transition is therefore first order. To this aim, we 
have performed a fit to the pseudo-critical couplings with the first order law 

(3 pc (L) =P C + ^, (31) 

finding 

(3 c (N t = A) = 3.245(3), B(N t = 4) = -500(30) , (x 2 /d.o.f. = 2.1) 

p e (N t = 8) = 5.852(8) , B(N t = 8) = -1300(100) , (x7 d -°- f - = 0.6) . (32) 

Looking at the x 2 /d.o.f., one can argue that for N t = 4 first order should be ruled out, 
whereas N t = 8 is compatible with first order scaling. H This can be due to the limited 
volumes [L < 150) considered for N t = 8 and to the larger error bars in the determinations 
of the P P cS with respect to the N t = 4 case. However, for N t = 8 the good agreement 
between the numerical result for the magnetic critical index and the corresponding value 
in the 2d XY model supports the claim that, even for this Nt, the transition is BKT. 

5 Conclusions and outlook 

The purpose of this paper has been to study the critical behavior of 3d U(l) LGT at 
finite temperatures, through the formulation on an asymmetric lattice. While the theory 
at zero-temperature is always in the confined phase, at finite temperatures it undergoes a 
deconfinement phase transition, just as it happens for 4d QCD. Analytical results from the 
high-temperature expansion suggest that this transition is of BKT type, but compelling 
numerical evidence is missing that indeed critical indices of 3d U(l) LGT coincide with 
those of the 2d XY model. 

This paper is the first step in the construction of the phase diagram of 3d U(l) LGT 
in the (j3 t , /3 s )-plane, where (3 S {(3 t ) is the spatial (temporal) coupling. In particular, we 
restricted ourselves to the case f3 s = and, by means of numerical Monte Carlo simulations 
on a dimensionally reduced effective theory, found evidence that the theory belongs indeed 
to the same universality class of the 2d XY model. The key observations have been the 
appearance of essential scaling and the agreement of the magnetic critical index rj with 
that from the 2d XY model. 

The next step is the extension of the numerical procedure established in this paper to 
the general case of (3 S ^ 0. 
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Figure 1: (Top) Absolute value of the Polyakov loop vs (3 on a 1 x L 2 lattice, with L=32, 
64, 128. (Bottom) Susceptibility of the Polyakov loop v s f3 on a 1 x L 2 lattice, with L=32, 
64, 128. For the L=128 case the multi-histogram interpolation around the peak is shown. 
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Figure 3: (Top) r/ cS for N t — 1 on lattices with several spatial sizes L at /3 = 1.12. For 
all lattices we fixed Rq=10. Errors are calculated with the jackknife method. (Bottom) 
L V C(R) versus R/L, with t] fixed at the central value of our determination through the 
method of the effective rj e s (see the text). 
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Figure 4: (Top) r] efr for _V t = 4 on lattices with L=64, 128, 200 at /3 = 3.42. For all lattices 
we fixed i?o=10. Errors are determined by the jackknife method. (Bottom) U I C{R) versus 
R/L, with r) fixed at the central value of our determination through the method of the 
effective r) e g (see the text). 
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Figure 5: (Top) r] cS for N t = 8 on lattices with L=64, 128, 200 at /3 = 6.38. For all lattices 
we fixed i?o=10. Errors are determined by the jackknife method. (Bottom) L V C(R) versus 
R/L, with r) fixed at the central value of our determination through the method of the 
effective rj e s (see the text). 



17 



